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ABSTRACT 



We present measurements of the shape of the stellar line-of-sight velocity 
distribution out to two effective radii along the major axes of the four elliptical 
galaxies NGC 2434, 2663, 3706, and 5018. The velocity dispersion profiles are 
flat or decline gently with radius. We compare the data to the predictions 
°f / = f{EiL z ) axisymmetric models with and without dark matter. Strong 
tangential anisotropy is ruled out at large radii. We conclude from our 
measurements that massive dark halos must be present in three of the four 
galaxies, while for the fourth galaxy (NGC 2663) the case is inconclusive. 



Subject headings: dark matter - galaxies: elliptical and lenticular, cD - galaxies: 
individual: NGC 2434, NGC 2663, NGC 3706, and NGC 5018 - galaxies: 
kinematics and dynamics - galaxies: structure 
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1. Introduction 

The presence of dark matter in normal elliptical galaxies is still controversial. Although 
none seems to be required inside the half-light radius R e , an extended halo of dark matter 
might surround the luminous regions (e.g., de Zeeuw & Franx 1991; Bertin & Stiavelli 
1993). Evidence for such dark halos has been found for some galaxies through studies of 
HI kinematics (e.g., Franx, van Gorkom & de Zeeuw 1994), X-ray emission (e.g., Jones 
et al. 1994), radial velocities of globular clusters (e.g., Mould et al. 1990) and planetary 
nebulae (e.g., Hui et al. 1994), or gravitational lensing (e.g., Maoz & Rix 1993). 

The shallow decline of the stellar line-of-sight velocity dispersion a seen in many 
normal ellipticals has also been taken as evidence for unseen mass at large radii (Saglia et 
al. 1993; Bertin et al. 1994; Carollo & Danziger 1994a, b, hereafter CDa,b). Unfortunately, 
a stellar orbital structure which becomes biased towards tangential motion at large radii 
can give rise to an almost flat a profile as well, and hence can mimic the presence of dark 
matter. This ambiguity in the interpretation of stellar kinematic data can be resolved by 
measuring the shape of the stellar velocity distribution along the line of sight (hereafter 
'velocity profile', VP) at large radii, since this depends strongly on the orbital anisotropy 
(Dejonghe 1987; Gerhard 1991, 1993). 

Accurate VP measurements in the outer, faint regions of galaxies are not easy to 
achieve, because a signal-to-noise ratio S/N > 50 per A is required to measure the lowest 
order anti-symmetric and symmetric deviations of the VP from a pure Gaussian (e.g., Rix & 
White 1992; van der Marel & Franx 1993, hereafter MF). Accurate VP shape measurements 
were obtained recently by, e.g., Bender, Saglia & Gerhard (1994) and van der Marel et 
al. (1994a, hereafter MRCFWZ), but these were restricted to radii < \R e - 

In this Letter we present VP shape measurements along the photometric major axes 
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of the four galaxies NGC 2434, 2663, 3706, and 5018, classified as E in the RC3 (de 
Vaucouleurs et al. 1991). They cover a range in luminosity, environment, metallicity and 
rotational properties. NGC 3706 and 5018 have faint shells (Malin & Carter 1983). The 
measurements extend to radii > 2R e , and allow us to quantify the anisotropy of the stellar 
velocity distribution and to determine whether dark halos are present. 

2. Velocity profile measurements 

Standard R ESO-Bessel broad-band imaging and deeply-exposed major axis long-slit 
spectra were taken for NGC 2434 and 3706 at the 2.2m telescope (EFOSC2; 4600-6000A; 
5A FWHM; exposure times 5 and 8 hrs., respectively) and for NGC 2663 and 5018 at the 
New Technology Telescope (EMMI; 4700-5500A; 4.5A FWHM; exposure times 6 and 5 hrs., 
respectively), at ESO, La Silla. Detailed information about these data and their reduction 
is given in CDa,b, who measured the i?-band surface brightness profiles and derived the 
mean line-of-sight velocity V and velocity dispersion a as a function of projected radius, 
based on the assumption of a Gaussian VP. 

We reanalyzed the spectra (rebinned spatially to have S/N > 50 per A at all radii) 
with the Fourier-Fitting package of Franx, Illingworth & Heckman (1989), as modified 
by MF to allow for non-Gaussian VP shapes. The VP is expanded into an orthogonal 
Gauss-Hermite series with parameters: (i) the line strength, mean V , and dispersion a 
of the best-fitting Gaussian VP; and (ii) the Gauss-Hermite moments /i 3 , /i 4 and higher 
orders that measure deviations from the best-fitting Gaussian. The parameter h 3 measures 
the lowest order anti-symmetric deviation from a Gaussian, /i 4 the lowest order symmetric 
deviation. A value /i 4 > generally indicates that the VP is more centrally peaked than a 
Gaussian, a value /i 4 < that it is more flat-topped. The S/N and instrumental resolution 
of our data are insufficient to measure higher order coefficients such as h 5 and h e . 
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The Fourier-Fitting program is based on an algorithm that minimizes the % 2 in the 
Fourier domain. For each galaxy, one template spectrum was constructed as a weighted mix 
of individual stellar spectra. The weights were chosen to optimize the fit to the observed 
galaxy spectra, using a technique similar to that of Rix & White (1992). Tests indicate 
that systematic errors in the VP parameters due to residual template mismatching and 
imperfect continuum subtraction are negligible for h 3 . For /i 4 they can be of order ±0.03, 
but they are generally independent of radius. 

The resulting radial profiles of V, a, /i 3 , and /i 4 for the four galaxies are shown in 
Figure 1. NGC 2434 has a rotating core, but shows negligible rotation beyond 5". The a 
profile has a central peak and is flat beyond 5". NGC 2663 has no or very small rotation, 
and has a steeply declining a profile. NGC 3706 and 5018 are fast rotators. NGC 3706 
has a central peak in cr, while NGC 5018 has a central minimum. For both galaxies, the a 
profile shows a modest decline in the outer parts. The measured velocity dispersions for all 
galaxies are significantly larger than the instrumental dispersion of ~ 130 km s -1 . 

The V and h 3 profiles generally have similar radial shape but opposite sign, in 
agreement with other studies (MF; Bender et al. 1994; MRCFWZ). In the outer parts 
of NGC 3706, where V/a is large, h 3 and V have the same sign. The VP asymmetries 
measured by h 3 are caused by the intrinsic skewness of the local stellar velocity distribution 
and by line-of-sight projection (e.g., Evans & de Zeeuw 1994). The /i 4 profiles differ 
considerably among the four galaxies. Out to > 2R e , /i 4 > for NGC 2434 and 5018, 
K £ - 0.05 for NGC 3706 and h A ^ - 0.1 for NGC 2663. 

3. Dynamical modeling 

We compare the data to axisymmetric dynamical models in which the phase-space 
distribution function (DF) / depends only on the two classical integrals of motion, the 
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energy E and angular momentum L z along the symmetry axis. The even part f e (E,L z ) = 
[f(E, L z ) + f(E, —L z )]/2 is determined uniquely by the luminous mass density Pl{R, z), for 
any potential ^{R, z) (R, (j) and z being the usual cylindrical coordinates). The odd part 
f Q can be freely specified, provided that / = / e + / Q > for all physical values of (E, L z ). 
Methods for constructing f(E,L z ) DFs for realistic mass distributions were developed only 
recently (Hunter & Qian 1993; Dehnen & Gerhard 1994; Qian et al. 1994, hereafter QZMH). 

The even and odd parts of the DF determine the even and odd parts of the VPs, respec- 
tively: VP e , (^) = [VP(u) ± VP(-v)]/2 = (l/£)J7J7e,o(£,£z) dv x , dv y , dz', where E is the 
projected surface mass density, (x',y') are on the plane of the sky, and z' is along the line of 
sight. Since we are interested in the mass distribution, we consider only VP e . We re-express 
the best-fitting Gauss-Hermite series as a Gauss-Hermite series with zero mean (van der 
Marel et al. 1994b). This yields: cr e , the dispersion of the best-fitting Gaussian to VP e , and 
Z4, the fourth Gauss-Hermite moment of VP e . The quantity a e is an observational estimate 
of crms, the RMS projected line-of-sight velocity. If VP e is Gaussian, then a e = crms- Fig- 
ure 2 displays the a e and z 4 profiles. If the VP is symmetric and there is no mean streaming, 
then a e = a and z 4 = /i 4 . This is the case for NGC 2434 and 2663 to within the observational 
errors. By contrast, in the fast rotators NGC 3706 and 5018 the a e profile declines more 
gently in the outer parts than the a profile, and the z 4 profile differs from the /i 4 profile. 

We consider models with pi, stratified on spheroids with axis ratio q, of the form: 

PL (R,z) = Po(jY(l + ^)", m 2 = R 2 + z>/q>. (1) 

The associated DFs and VPs were discussed by QZMH. The isophotes of the projected 
mass density have ellipticity e = 1 — q', where q' 2 = cos 2 i + q 2 sin 2 i, and i is the inclination. 
For each galaxy, p , b, a, and (3 were chosen to fit the observed i?-band surface brightness 
profiles, after convolution of the model with the seeing. Since we are interested in the outer 
regions, e was fixed to the average observed value outside 0.5R e . It is listed in Table 1, 
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together with the effective radii and the best-fitting parameters b, a, and (3. These are 
independent of i. 

Since the calculation of f(E,L z ) from pi, requires substantial numerical effort, we 
adopted two simplifications. First, we use the Jeans equations to calculate the second 
velocity moments (vjj) and Since (vjj) = {v 2 z ) and (vrv z ) = in the two-integral case, 

crms follows by line-of-sight integration at projected radius R' . We compare the resulting 
crms to the observed a e . The errors thus introduced are small, because the observed z 4 
values are relatively close to zero. Second, we only calculate z 4 in the limit of very large 
radii (R/b ^> 1), where its calculation is much easier (QZMH). The observed z 4 values 
extend to R'/b % 5 in NGC 2434 and 2663, and to R'/b % 2 in NGC 3706 and 5018, where 
R' is the projected radius on the sky. We verified that these simplifications do not influence 
any of our main conclusions, by calculating the full DF for a number of specific cases. In 
particular, the z 4 values for the range R'/b % 2 — 5 are very close to the large radii result. 

Figure 2 shows the predicted crms profiles for the self-consistent case (i.e., no dark 
halo), for edge-on models (i = 90°). The normalization is determined by the average 
mass-to-light ratio of the stellar population M/L, which can be chosen freely. Independent 
of the assumed M/L, we conclude the following. For NGC 2434 the predicted profile falls 
more steeply in the outer parts than the observed profile, while for NGC 2663 it falls less 
steeply. For NGC 3706 and 5018 the models predict too little motion in the outer parts 
relative to the inner parts. These results agree with those of CDa,b, who discussed also the 
effects of varying i and of adding a dark halo. 

The dashed curves in Figure 3 show the predicted z 4 values in the large radii limit as 
function of i. The dependence on i is weak, except at the smallest allowed inclinations. 
As expected, the predicted z 4 is always strongly negative, because any flattened f(E,L z ) 
model is supported mainly by azimuthal motion (e.g., Dehnen & Gerhard 1994). The 
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hatched bands represent the error bars on z 4 at the outermost measured point, and are a 
conservative estimate of the observational uncertainty in the mean z 4 values measured at 
R' > R e . Except for NGC 2663, the predicted z 4 values lie well below the observed values. 
This is not caused by the simple form adopted for pi, and remains so when a, /3, and b are 
varied within the range allowed by the errors of the best fit (Table 1). It also cannot be due 
to systematic errors in the observations, as these are |Az 4 | <^ 0.03 (of the same order as for 

K). 

Models with f(E,L z ) and a dark halo can fit the observed a e profiles reasonably 
well (CDa,b). We study the z 4 values predicted by f(E,L z ) models with a dark halo, by 
(non-selfconsistently) embedding the luminous mass-density pi in a power-law potential 
^dk oc (R 2 + £ 2 /<?dk 2 ) 7 ^ 2 at large radii, which corresponds to a total density profile in the 
equatorial plane p oc i? 7-2 (Evans 1994). We use again the asymptotic results of QZMH for 
R/b ^> 1. In this limit a e oc R n ^ 2 . For each galaxy, 7 was fixed by fitting to the observed a e 
profile at large radii. The solid curves in Figure 3 show the predicted z 4 for different values 
of the axis ratio gak (= 0.8,0.9, 1.0) of the potential ^dk- These correspond to axis ratios 
9mdk of the halo mass density ranging between 0.4 and 1.0. The predicted z 4 increases when 
<7dk decreases, since fewer high angular momentum orbits are needed to reproduce the same 
flattened mass density in a flatter potential. For NGC 2434 the z 4 values overlap with the 
observed values for q m dk ~ 0.4 — 0.8. However, for NGC 3706 and 5018 very flattened halos 
(<7mdk ~ 0.1 — 0.4) are required to bring the predicted z 4 into agreement with the observed 
value. 

4. Discussion and Conclusions 



The VPs of elliptical galaxies can now be measured reliably out to ~ 2R e . We have 
presented the first such measurements for four galaxies, and have compared them to the 
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predictions of axisymmetric f(E,L z ) dynamical models. Radially anisotropic models 
generally predict VPs that are more centrally peaked than a Gaussian, while tangentially 
anisotropic models predict VPs that are more flat-topped (e.g., Gerhard 1993). The 
observed z 4 values imply a velocity distribution that is less tangentially anisotropic than in 
a self-consistent f(E,L z ) model. However, such a velocity distribution predicts a a e profile 
that falls off steeper than in the f(E,L z ) model (e.g., Dehnen & Gerhard 1993). This is 
consistent with the data for NGC 2663, but it increases the discrepancy with the gently 
declining a e profiles of NGC 2434, 3706, and 5018. Thus we conclude that for NGC 2434, 
3706 and 5018, no self-consistent axisymmetric model can fit simultaneously the a e and the 
VP data. These galaxies must have a dark halo. 

We have shown that f(E,L z ) models with flattened massive dark halos predict 
z 4 -values in harmony with the measurements for NGC 2434, 3706 and 5018. The f(E,L z ) 
dark halos for the latter two galaxies are very flat. Although we cannot exclude this 
possibility, we suspect that NGC 3706 and 5018 are surrounded by less-flattened halos. 
The resulting discrepancy in z 4 can then be explained either by small systematic errors 
in the z 4 measurements, or by an orbital structure which is less tangentially anisotropic 
than our f(E,L z ) models. So, while for NGC 2434 possibly / = f(E,L z ), for NGC 3706 
and 5018 it seems more plausible that / ^ f(E,L z ). To give these qualitative arguments 
a more quantitative basis, and to determine the size and mass of the dark halos, it will be 
necessary to construct flattened three-integral models for the galaxies in our sample. 

We thank Martin Schwarzschild and Wyn Evans for comments on an earlier version of 
this paper. During part of this work, EEQ was supported by NSF grant DMS-9407559, 
and RPvdM was supported by NASA grant HF-1065.01-94A, awarded by the STScI which 
is operated by AURA, Inc., for NASA under contract NAS5-26555. 
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Table 1: Basic Properties and Model Parameters' 1 



Galaxy 


Re(") 


e 


a 


P 


K") 


NGC 2434 


24 


0.08 


-1.9 


-0.52 


13 


NGC 2663 


50 


0.30 


-1.6 


-0.55 


10 


NGC 3706 


27 


0.35 


-2.3 


-0.53 


31 


NGC 5018 


22 


0.30 


-2.1 


-0.90 


31 



"Listed are: the effective radius R e (Lauberts & Valentijn 1989), the average ellipticity e 
outside R e /2 (CDa,b), and the model parameters a, /3, and b of equation (1). Mean errors 
on the best fit parameters are Aa = ±0.1, A/3 = ±0.05, and Ab = ±5". 
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Fig. 1. — Velocity profile (VP) parameters for the four galaxies in our sample as functions of 
projected radius R' on the plane of the sky. From top to bottom: the mean V and dispersion 
a of the best-fitting Gaussian to the VP, and the Gauss-Hermite coefficients h 3 and /i 4 . 
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Fig. 2. — The dispersion a e of the best-fitting Gaussian to the even part of the velocity 
profile, VP e , and the Gauss-Hermite moment z 4 of VP e , as functions of projected radius 
R' for the four galaxies in our sample. The curves in the panels for a e are the (seeing 
convolved) predictions for ormSj the RMS projected line-of-sight velocity, for edge-on models 
with distribution function f(E,L z ) and no dark halo, as described in the text. 
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Fig. 3. — Behavior of z 4 as a function of inclination angle i for our (a,/3) models for the 
four galaxies in our sample. A value z 4 < indicates that VP e , the even part of the velocity 
profile, is more flat-topped than a Gaussian; a value z 4 > that it is more centrally peaked. 
Dashed curves: f(E,L z ) models without dark halo; Solid curves: f(E,L z ) models with 
a power-law dark halo with, from top to bottom, axis ratio of the isopotential spheroids 
9dk = 0.8, 0.9, and 1.0. The hatched regions indicate where the z 4 data points with R' > R e 
lie. 



